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1 Abstract 

In this paper we discuss the dynamical features of intermittent fluctuations in homo- 
geneous shear flow turbulence. In this flow the energy cascade is strongly modified by 
the production of turbulent kinetic energy related to the presence of vortical structures 
induced by the shear. By using direct numerical simulations, we show that the refined 
Kolmogorov similarity is broken and a new form of similarity is observed, in agreement 
to previous results obtained in turbulent boundary layers. As a consequence, the in- 
termittency of velocity fluctuations increases with respect to homogeneous and isotropic 
turbulence. We find here that the statistical properties of the energy dissipation are prac- 
tically unchanged with respect to homogeneous isotropic conditions, while the increased 
intermittency is entirely captured in terms of the new similarity law. 
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2 Introduction 



The issue of anomalous scaling laws in turbulence has been largely addressed for the 
idealized case of homogeneous and isotropic turbulence once it became clear that a purely 
dimensional power law [jlj can not be consistent with the intermittent nature of energy 
dissipation Pf. The Kolmogorov-Obukhov refined similarity (RKSH) P] provides the 
key theoretical point by relating the anomalous correction of the scaling exponents of 
velocity structure functions directly to the statistical properties of the dissipation field. 
The RKSH has been subject to close scrutiny by a number of independent investigations 
making use of both experiments and direct numerical simulations and it may be considered 
a well assessed physical law. One of the early difficulty for an accurate evaluation of 
scaling exponents, associated to the existence of too small scaling ranges in moderate 
Reynolds number turbulence, has been recently overcome with the use of the Extended 
Self-Similarity (ESS) [|J. This approach has allowed accurate estimates of the exponents 
from moderate Reynolds number flows || and the evaluation of the RKSH, in its extended 
form, from DNS data 

A more recent research interest in turbulence has been focused on the assessment of 
RKSH in non-isotropic and non-homogeneous conditions, namely in flows characterized 
by strong mean shear and anisotropy. The main focus is to understand to what extent 
RKSH may be considered still able to capture, in such conditions, the essential features 
of turbulence dynamics. 

In this context recent results for a turbulent channel flow |7]], have pointed out how 
the classical RKSH holds in the higher part of the logarithmic region. This result should 
not be surprising since the RKSH is established as a consequence of the balance between 
energy transfer due to non linear interactions and dissipation. In the log-layer, where 
the mean shear is weak, the dynamics of the flow is characterized by the inertial energy 
transfer. In fact, when approaching the wall, the shear becomes larger and larger and 
a clear failure of the RKSH is observed J7|. This result gives a partial answer to the 
point raised above since it gives evidence that, in the near wall region, the homogeneous 
isotropic behavior is far from being recovered. 

At this point a question naturally arises about the existence of an alternative scaling 
law in such conditions. On the basis of physical considerations on the respective role of 
energy production and energy transfer a different Kolmogorov-like scaling law has been 
proposed in [0] and preliminarily checked against DNS data from the buffer region of 
a turbulent channel flow. The validity of this similarity law has been confirmed by a 
successive experimental work in a turbulent boundary layer. However, certain aspects 
of the problem are still not completely clear. Actually, the new similarity law has been 
conceived as a direct consequence of the presence of a strong shear, but it has been verified 
only in wall turbulence where, in principle, some other features such as the continuous 
variation of the local shear, the non uniform momentum flux across adjacent layers and 
the suppression of wall normal velocity fluctuations may play a quite significant role. 
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To isolate the effect of the shear we reconsider here the problem in the context of a DNS 
of a homogeneous shear flow in a confined box. In fact, in the rationale of dimensional 
analysis, this flow presents the minimum level of complexity though maintaining the 
essential features we want to address: thus, it is ideal to study the effect of a pure shear 
avoiding other concurrent effects. On the other hand, it allows to exploit homogeneity 
in all spatial directions so that a more complete and accurate statistical analysis can be 
performed, at a level of detail that could never be achieved in the very limited buffer 
region of a wall bounded flow. 

As a counterpart, the homogeneous shear flow presents certain drawbacks related to 
the artificial nature of its confinement. To understand to what extent the present flow may 
serve as a prototype for shear-dominated flows, we do analyze in detail the dynamics of the 
coherent vortical structures observed in the numerical simulation. This is instrumental to 
qualify the flow, i.e. its dependence on the aspect ratio of the computational box, in view 
of the use we do of the flow itself for the evaluation of scaling laws and their assessment. 

Homogeneous shear flows have been investigated in the literature in many theoretical, 
numerical and experimental works. Rogers and Moin || and Lee et. al. |J studied the 
topology and mutual interactions of vorticity showing that many dynamical features are 
quite similar to those observed in the wall region of a turbulent boundary layer. A more 
detailed study on the same subject has been performed by Kida and Tanaka [|10| who 
discussed the regeneration cycle of the streamwise vortices in a homogeneous shear flow. 

Experimentally the papers by Rose jy]], Champagne Harris and Corrsin [12|, Tavoularis 
and Corrsin [TJj], have investigated homogenous shear flows obtained in a wind tun- 
nel, focusing on the local isotropy of small scale fluctuations. Following the Kolmogorov 
theory, an important question to investigate, which becomes extremely well posed in the 
case of the homogeneous shear flow, concerns the recovery of small scales local isotropy in 
turbulent flows characterized by large scale non isotropic forcing. The same kind of ques- 
tion has been addressed by Saddoughi and Veeravalli |15[] who investigated local isotropy 
for the turbulent flow in a logarithmic boundary layer. 

The issue of local isotropy has been also addressed more recently by Pumir and 
Shraiman [|1^] and by Pumir |17]] who performed an extensive analysis of numerical simu- 
lations of statistically stationary homogeneous shear flows. Interestingly enough, most of 
the statistical properties found by Pumir [17] are quite close to the experimental findings 



of Garg and War haft [|18|] and Shen and Warhaft fll9fl . In the latter experiment, for the 
first time, an active grid was used in order to produce a homogeneous shear flow in a wind 
tunnel whose integral scale in the streamwise direction is rather constant. 

Taking full advantage of most of these results, and in order to be able to reconsider 
the issue of local isotropy from a different perspective, as we plan to do in the near 
future, our interest here is mainly focused on the relationship between intermittency and 
the anisotropy induced by the mean velocity gradient at large scale. In the context of 
boundary layer we have shown how the increase of intermittency may be explained in terms 
of the phenomenology described by the new form of scaling law. In particular, as suggested 
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in 0, a new length scale L s = yt/S 3 enters in the description of the statistical properties 
of turbulence, where e is the mean rate of energy dissipation and S the mean large scale 
shear. For scales smaller than L s the statistical properties of the turbulent fluctuations 
are similar to those observed in homogenous and isotropic flows while for scales larger 
than L s the refined Kolmogorov similarity (RKSH) is broken and intermittency increases. 
The phenomenological theory developed for the turbulent boundary layers relies upon 
the mean shear strength S as the only parameter which fixes the scale L s where a new 
form of RKSH should be observed. For this reason, the homogenous shear flow becomes 
a natural test case to understand whether or not the phenomenological theory devised 
for the boundary layer is sufficiently general to be applied to any turbulent flow. As we 
shall discuss in this paper, our numerical simulations of homogenous shear flows support 
rather well the new phenomenology, leading to possible important implications in our 
understanding of shear turbulence. 

The paper is organized in the following way. In section 3 we briefly describe how the 
numerical simulations are performed. In section 4 we discuss the regeneration cycle of the 
vortical structures observed in the flow while in section 5 we address the intermittency 
cycle which characterizes the dynamical behavior of the system. In section 6, after a 
short review of the phenomenological theory proposed for the turbulent boundary layer, we 
describe the intermittency properties of turbulent fluctuations. In section 7, we investigate 
the anomalous scaling of the energy dissipation field, and, finally, in section 8 we draw 
our main conclusions. 

3 Homogeneous shear flow 

We consider here a turbulent flow in a confined box with an imposed mean shear S free 
from boundaries. The full Navier-Stokes equations are solved (DNS), after decomposing 
the velocity field into mean value and fluctuation 



where x G V = [— X X ,X X ] X [— X y , X y ] x [— X Z ,X Z ], V identifies the computational box, 
U {%/) = S(y + X y ) is the mean flow and e\ is the unit vector in the streamwise direction x. 
The mean gradient S is in the normal direction y while z denotes the spanwise coordinate. 
The Navier Stokes equations are written in terms of velocity fluctuations 



where ( is the vorticity, v the normal velocity, p the pressure and v the kinematic viscosity. 
To achieve efficiency and accuracy a Fourier spectral method is advisable to solve the 



v(x, t) = U (y)efi + u(x, t) , 



(1) 




du 
dx 



(2) 
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initial value problem for eq. @. However, due to the mean flow, the intrinsically non- 
periodic term U(y)du/dx appears in the equations. Luckily, by using Rogallo's technique 
||20||, this difficulty is removed by the transformation 



(3) 



In the new variables periodic boundary conditions can be enforced in all spatial directions 
allowing efficient pseudo-spectral methods to be used for spatial discretization. The time 
integration is performed using a third order low storage Runge-Kutta method |2]J and the 




nonlinear terms are fully de-aliased by zero padding. Since the image of the computational 
box in physical space gets distorted, eq. ([|), a re- meshing procedure is periodically applied 
to allow long time integrations. Using periodicity in the £ direction, the computational 
domain is transformed back into a non skewed domain whenever the plane y — X y has 
moved by 2X X in the streamwise direction, i.e. every At r = X x /SX y . This procedure 
may introduce aliasing errors since the dynamics of wave vectors in physical space is 
time dependent [p^j , as seen from the relation between the Fourier transforms in the 



computational and in the physical space, 

Ui{k x , k y , k z , t) = Ui{k^ k v - Srk^, k r). (4) 

To avoid aliasing errors, the spectral components with wave numbers outside the interval 
[—kmax, k max ] are set to zero at re- meshing. This filtering introduces a characteristic wave 
number defined as 

h 

7 n max z_\ 

k r = . 5 

where k max = 7rN y /2X y and A = X x /X y is the aspect ratio of the computational domain. 
In order to asses the influence of the filtering, the aspect ratio of the box and the resolution 
of the grid were varied. In all cases k r was larger than the dissipation wavenumber and 
no appreciable alteration of the dynamics was observed. 

A crucial point for our successive analysis is the ability of the system to reach statistical 
stationarity. Under this respect, our computations entirely confirm the results of Pumir 
and Shraiman |16[ and Pumir [jT7|j . Actually the evolution of the turbulent kinetic energy 



is described by the equation 

|[y] + ^H = -WC 2 ] (6) 

where the square brackets denote spatial average. This equation suggests a possible 
balance between the production S[uv] and the dissipation v[(, 2 }- This balance is not 
achieved for the instantaneous fields, as shown in figure 1 which reports the history of 
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turbulent kinetic energy in our longest calculation. Large fluctuations in energy (42%) 
and enstrophy (50%) are apparent in the pseudo- cyclic behavior shown in the figure. 
The fluctuation level in this system is quite huge, expecially when compared to that of 
forced homogeneous isotropic turbulence, which is limited to only a few percent of the 
rms value (5 — 6%) for both kinetic energy and enstrophy. In fact, the pseudo-cyclic 
behavior corresponds to statistically stationary conditions, as follows from time averages 
performed over time periods much larger than the typical length of the cycle. 

Clearly, DNS enables to achieve a good level of homogeneity, as shown in figure 2. In 
particular both the turbulent kinetic energy and the Reynolds stresses are substantially 
constant in the y direction. 

In order to check the dynamics of the velocity gradients, the probability density func- 
tion of du/dx and of du/dy have been computed, figure 3 and 4. They are often used to 
asses the small scale dynamics of the flow and to characterize its degree of anisotropy. The 
computed values of skewness and flatness are in good agreement with the experimental 
results at Re\ ~ 100 of Shen and Warhaft |19| and of Ferchichi and Tavoularis [|3[] at 
i?e A = 140. 



4 Regeneration cycle of vortical structures 

The energy fluctuations discussed in the previous section correspond to a regeneration 
cycle of the vortical structures. Most of the previous numerical simulations have described 
the main features of coherent structures and the vorticity statistics for the early stages 
of development of the flow. In particular, Rogers & Moin || showed typical hairpin 
vortex structures remarkably similar to those observed in wall bounded flows. Lee el. 
all. |§ discussed streamwise vortices and high and low speed streaks commenting on the 
similarities with the buffer region of wall bounded flows. In a more recent work Kida & 
Tanaka [l(| proposed a regeneration mechanism for the streamwise vortices and pointed 



out the role of vortex sheet instability in the formation of new vortices. The present 
analysis deals with the same issues, but it is focused on the statistical steady state regime 
of the flow. 

The phenomenology is better understood by considering a few snapshots along the 
history of the turbulent kinetic energy and of the Reynolds stresses. In figure 5 the 
large energy bursts already noticed in the previous section are clearly correlated to large 
negative values of the Reynolds stresses. The bursts are induced by large rates of energy 
injection from the mean flow associated with the presence of streamwise vortices. Figure 
6, corresponding to an instant at the beginning of a stage of energy growth, actually shows 
a large population of quasi-streamwise vortices, here visualized through the discriminant 



of the velocity gradient ||24j| . During the successive phase of energy growth the streamwise 
vortices, by interacting with the mean velocity field, give rise to instantaneous profiles 
characterized by the typical ramp and cliff pattern. The pdf of du/dy is a suitable tool 
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to characterize statistically the ramp and cliff structures, figure 4. The skewness of 0.82 
gives reason of intense positive events of du/dy much more probable than negative ones. 
As for the passive scalar, the streamwise velocity recovers the mean gradient through 
ramps followed by cliffs ||25| . The ramps correspond, on average, to streamwise velocity 
increasing with y, hence a resulting negative spanwise vorticity is expected, on average. 
The pdf of the spanwise vorticity, figure 7, confirms this analysis. 

Ramp and cliffs are associated with thin regions of (negative) spanwise vorticity, i.e. 
vortex sheets, as shown in figure 8. The sheets become unstable and eventually roll up into 
spanwise vortices through the classical Kelvin-Helmholtz mechanism, see figure 9 where 
spanwise vortices are systematically located in the regions where the roll-up process is 
occurring. Vortex sheets and spanwise vortices are clearly responsible for the observed 
energy bursts by producing large negative Reynolds stress events. 

Immediately after each burst, the non- linear interactions are enhanced and the original 
vortex structures lose their order resulting into a randomized vorticity field. This phase, 
shown in figure 10, reduces the level of anisotropy of the flow. However, soon afterwards, 
the shear term of the vorticity equation enforces a mean orientation to the structures, see 
figure 11, corresponding to a minimum for the energy, just before the occurrence of the 
successive burst. The vortical structures are now getting more and more aligned with the 
streamwise direction and a new cycle starts. 

The regeneration cycle just analyzed here for the statistical steady state phase is 
very similar to that described by Kida & Tanaka |I| for the early stages of evolution of 
the flow from isotropic initial conditions, suggesting that the operating mechanisms are 
substantially identical. 



5 Intermittency cycle 



The increase of the turbulent kinetic energy due to the streamwise vortices can be under- 
stood, basically, in terms of the linear lift-up mechanism and the related transient growth. 
A non linear mechanism is required instead to explain the saturation and the break-down 



of the ordered system of vortices. In spectral space, as already discussed by Pumir [T? 
the Fourier mode (0, 0, ±1) gives the leading contribution to the energy growth. From the 
linearized equation for this mode, 



du 
~dt 

dv 
~dt 



-Sv — vu 



—isv : 



(7) 



a growing amplitude is expected whenever S[Re{u)Re{v) + Im(u)Im{y)} < 0. As soon 
as the energy grows appreciably, the energy drained by the interactions with the other 
modes in favor of the smaller scales quite rapidly saturates its amplitude. Hence, the 
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characteristic frequency of the intermittency cycle is determined by the dynamical balance 
between growth of the basic mode and energy transfer. 

In principle, if the mode (0, 0, ±1) were isolated from the others during its growth, 
the energy due to the interaction with the mean shear would be almost entirely found in 
this mode. Then, the characteristic time of the system would essentially coincide with 
the eddy turn over time of the basic mode, hence to the spanwise size of the box L z . 

In fact, the leading mode evolves in presence of many others. In this case the saturation 
time, i.e. the time of the effective activation of the nonlinear energy transfer, is related 
to the separation, in wavenumber space, between the leading mode and the small scales 
modes. The activation time becomes shorter and shorter as the distribution of energy 
among all Fourier modes approaches a continuous distribution. For the statistical steady 
state regime, in particular, the non linear energy transfer is crucial from the very beginning 
of the growth. Actually the energy spectrum in correspondence of the large bursting 
phenomena, figure 12, shows initialy a growth which, though dominated by the mode 
(0,0, ±1), is spread all over the modes in the first decade. The evolution of the spectra 
manifest a strong energy transfer which though occurring during the entire cycle, is clearly 
prevailing in the phase of energy decrease. Hence the bursting frequency can not be 
estimated from the isolated dynamics of the sole (0, 0, ±1) mode, and, as a consequence, 
there is no a priori reason for a strong dependence of the bursting period on the size of 
the box. 

In view of a quantitative analysis of this issue, let Ejj and E u denote the kinetic energy 
of the basic flow and of the fluctuations, respectively. Clearly, E u is a function of time. 
Two quite different cases can be distinguished, namely 

a) E u > E v at t = 

b) E u < E v at t = . 

Case b) corresponds to a problem of transition to turbulence from small, though finite, 
disturbances and will not be considered here in further detail. We are presently interested 
in case a), which is more appropriate to describe statistical stationarity. Here t — should 
be understood as an arbitrary instant of time along the cyclic evolution of the system. 

Given the correlation function C(r) =< E u (t + r)E u (t) >, let us consider the correla- 
tion time, Tc, defined, e.g., as the smallest r such that C(tc) = 1/2C(0). Following the 
discussion of section 4, tq corresponds, roughly, to the characteristic time of production 
of large scale velocity/vorticity instability which is the primary forcing mechanism of the 
homogeneous shear flow. If r c is almost independent or, at least, weakly dependent on 
L z , the dynamical behavior of the system may be considered as substantially unaffected 
by the finite size of the box. A rigorous analysis of the correlation time dependence on 
the spanwise size L z would require a large number of highly resolved numerical simula- 
tions far beyond the present computer capabilities. In figure 13, we show a preliminary 
evidence that the behavior of the kinetic energy is not too sensitive to L z , by comparing 
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the simulations performed with L z = 2ir and L z = Air. The characteristic time of the 
generation mechanism of large scale vorticity seems not to depend on L z . These results 
suggest that the finite size of the box does not affect in a significant way the dynamics 
and the statistical properties of the homogeneous shear flow, despite the importance of 
the basic mode, whose length scale coincides with that of the box, in the dynamics of the 
flow. 



6 Intermittency and scaling 

In the previous sections we have analyzed the influence of the shear on the dynamics of 
the vortical structures. We are now ready to discuss its effect on statistical features of 
turbulent fluctuations, such as scaling laws and intermittency. 

Concerning the velocity field, a quantitative description can be given in terms of the 
possible scaling behavior of the structure functions, i.e. the moments of longitudinal 
velocity increments 

< SV P >=< {[u(x + r,t)-u(x,t)] ■ -} p > (8) 

where angular brackets denote ensemble averaging. For homogeneous and isotropic con- 
ditions, the dimensional prediction of Kolmogorov theory (K41) provides a scaling law in 
terms of separation Jl| 

<SV P > oc e p/3 r p/3 (9) 

where e denotes the mean rate of energy dissipation per unit mass. For p = 3, equation 
(|) is consistent with the exact result usually referred to as the four fifth law, 

<5V 3 > = —er. (10) 
5 

This equation has been widely confirmed by a number of experimental investigations ||] 
which, at the same time, showed that dimensional scaling does not hold for moments 
different from the third. In fact, a more complex dependence of the scaling exponents is 
found, 

<5V P > oc r cp , (11) 

with £ p a nonlinear convex function of p. This anomalous scaling related to intermittency, 
has been actively investigated over the last twenty years. 
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6.1 Dimensionless parameters 

To understand how intermittent fluctuations are affected by the shear, it is worthwhile to 
consider the typical length scales involved in the process and the relevance of the shear 
production term over the nonlinear inertial interactions. In presence of shear, velocity 
fluctuations arise for two main reasons, advection of streamwise momentum across the 
mean gradient and non-linear mixing. Let us denote by 5U = Sr the velocity difference at 
scale r due to the mean flow and by Su the fluctuation at the same scale. A rough estimate 
for the latter is Su oc e^V 1 / 3 . The length scale for which turbulent fluctuations are of 
the same order of magnitude as the velocity increments induced by the shear, uniquely 
identifies the shear scale L s 

£,= J|, (12) 



s 



that separates two different sub-ranges within the inertial range. For L ^> r ^> L 
(with L the integral scale), SU 3> Su and the dynamics and statistics of turbulence is 
expected to be dominated by the momentum flux due to the Reynolds stresses, i.e. by the 
production of kinetic energy. For rj <C r <C L s , Su ^> SU and the fluctuations induced by 
the mean shear are negligible with respect to those typically produced by the nonlinear 
term. In this case the dynamics of the flow is characterized by the energy transfer due to 
the same inertial interactions which are typical of isotropic turbulence. 

In order to quantify the predominance of either mechanism over the other, we consider 
the dimensionless parameter, defined as the ratio of an inertial and a shear time scale, 

ST? 2 

S* = ^j- (13) 
with q 2 =< UiUi >. S* can be interpreted as the ratio of two length scales 

/ , \ 2/3 

sr = (£) (14) 

where Id = g 3 /e. It follows that for large values of S* the dynamics of most of the inertial 
scales is dominated by the mean shear. 

The other dimensionless parameter needed for a complete description of a homoge- 
neous turbulent shear flow can be given as the ratio of the dissipative and the shear time 
scale, 

S: = S{v/efl\ (15) 
or, equivalently, as the ratio of the Kolmogorov length and the shear scale 
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Since S* measures the separation between the shear scale and the dissipative range, a 
sub-range of isotropic behavior may be recovered within the inertial range when it is 
small. 

Typical values of S* and S* are shown in figure 14, which is a compilation of several 
independent numerical and experimental simulations. The variety of behaviors observed 
in these simulations can easily be interpreted in terms of the location of the shear scale 
L s , quite different from one case to the other. 

6.2 Intermittency 

Let us now turn to the issue of intermittency. A quantitative measure of intermittency 
is provided by the flatness factor of the random variable 8V(r). Figures 15 and 16 show 
the probability distribution function of 5V(r) for two different separations both in the 
shear flow and in homogeneous isotropic turbulence. For small separations in both cases 
the pdf is clearly non Gaussian and the flatness factor, shown in figure 17, grows as 
separation decreases confirming that the dimensional scaling of the structure functions 
is clearly violated. However, important differences are apparent. The tails for the shear 
flow are quite higher and the pdf is more skewed towards negative values. For the larger 
separation, instead, the two pdfs are much more similar, though still not Gaussian. The 
Gaussian behavior is recovered only for r comparable with the integral scale. Globally, 
we find a correspondence with what was found in the turbulent channel flow j7|, where at 
small separation the pdfs of 5V(r) in the buffer and in the log-region differs substantially, 
with the negative tail much higher than the positive and an increased flatness factor in 
the buffer. 

6.3 Similarity laws 

The issue of intermittency can be quantified by looking directly at the scaling properties of 
the structure functions. As well known, starting from Landau objection || a revised form 
of similarity law has been proposed by Kolmogorov and Obhukov 0] (K62): assuming 
that the dissipation field is a random variable, equation (^) is corrected to 



where < e q . > denotes the q-th moment of the local energy dissipation field e\ oc averaged 
over a volume of characteristic dimension r. Given for the dissipation field a scaling law 
in terms of separation, 



< 5V P > oc < eP /3 > r 



.P/3 



(17) 



< e q r > oc r 



(18) 



from equation ( |ITD the scaling exponents of the structure functions are expressed as 



a P ) = r(p/3)+p/3 



(19) 
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with r(p/3) a nonlinear function of p and r(l) = 0. Equation fll9|) is consistent with 
the intermittent nature of the random variable 5V(r), i.e. its increasing flatness with 
decreasing separation, since the flatness factor, within the inertial range, turns out to be 
an unbounded function, diverging for small separations as r^ 4-2 ^ 2 . The intermittency is 
related to bursty signals characterized by localized and intense fluctuations and to the 
presence of coherent structures and associated regions of intense gradients which give rise 
to a highly non uniform dissipation field. 

Scaling laws in terms of separation are found only for sufficiently large values of the 
Reynolds number i.e., so far, only in experimental facilities. 

An extension of the scaling range can be achieved by using the Extended Self Similarity 
(ESS) as proposed by Benzi et. al. M, who introduced a generalized form of the scaling 
law(O), 

< 5V P > oc < 5V 3 >^ /C3 . (20) 

Using experimental data, the scaling law ( p0| ) has been shown to hold even at low and 
moderate Reynolds number, i.e. in cases, such as those amenable to DNS, where the scal- 
ing ([□]) is not directly detectable. Most interestingly, the ratios Cp/C3 have systematically 
been found independent of the Reynolds number and consistent with the measurements 
in the high Reynolds number regime. The ESS allows a generalization of the refined 
Kolmogorov similarity hypothesis ||, eq. (|T7|) , namely 

< e p l 3 > 

< 5V P > oc r r, < 5V 3 > p/3 . (21) 

This equation, which we will hereafter refer to as RKSH, has been carefully checked 
against experimental data for homogeneous and isotropic turbulence and confirmed also 
in those cases where neither < e v J 3 > nor the structure functions showed a scaling behavior 
in r. The two terms in eq. (|2"T|), namely < e p J 3 > and < 5V 3 >, take into account the 
fluctuations in the energy dissipation rate and the energy transfer through the inertial 
scales, respectively. 

Scaling laws as discussed so far heavily rely on the assumption of isotropy. In many 
turbulent flows, however, isotropy is broken by the mean shear and by the boundary 
conditions and a large production of turbulent kinetic energy occurs due to the interaction 
of the mean flow with the turbulent fluctuations. We may wonder if scaling laws emerge 
also in such conditions, and, in case, under which respects they differ from the classical 
ones. 

Recent numerical investigations of wall turbulence have shown how intermittency in- 
creases approaching the wall and how the scaling exponents ( p are not consistent with 
equation (|2"T|). Actually, in the wall region, the main dynamical process is represented by 
the momentum transfer occurring in the wall normal direction associated with a large pro- 
duction of turbulent kinetic energy via the Reynolds stresses dll/dy < uv >. In this case 
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a term of the form dU/dy < 5V 2 > is expected to account for the new balance between 
production and dissipation in the Karman-Howarth equation. Based on this analysis a 
new form of similarity law has been proposed for the wall region in terms of the structure 
function of order two 

< e p / 2 > 

<5V P > oc r M < 5V 2 > p/2 . (22) 

As discussed in @ for a turbulent channel flow, the classical RKSH, eq. (plf ), holds 
at the center of the channel in the bulk region where the momentum transfer due to the 
shear term is negligible with respect to inertial transfer. The new form of similarity law 
holds in the near wall region where high turbulent kinetic energy production occurs. 

We are presently going to check the consistency of equation (p2|) in the homogeneous 
shear flow where, in the absence of solid walls, the mean shear S is constant. In this flow 
the effect of the shear is isolated from other concurrent effects such as the suppression 
of wall normal velocity fluctuations and the non uniform momentum transfer across ad- 
jacent layers in the normal direction. Moreover the statistical analysis is simplified since 
homogeneity is exploited in all spatial directions. 

In figure 18 the similarity law for p = 6, is applied to DNS data for both isotropic 
and homogeneous shear turbulence. It holds for all the resolved scales in the isotropic 
case. For the shear flow, a clear failure is observed instead. In this case two different 
slopes may be extracted from the data, one in the dissipative region with the trivial value 
of 0.99, the other, in the inertial region, with value of 0.92. Analogous results are found 
for p = 9, figure 19, with a single slope for isotropic turbulence. For the shear flow, the 
dissipative slope is 0.98 confirming the quality of our DNS while the inertial fit gives a 
slope of 0.88. 

A similar analysis is presented in figures 20 and 21 for the new similarity law, eq. 
(f22j), in the case of the homogeneous shear flow, for p = 4 and p = 6. A unique slope 
can be effectively fitted all over the range of resolved scales with the same value of 1 
independently of the order of the moment. The failure of the classical RKSH and the 
validity of the new form for the shear flow are still better appreciated in the compensated 
plot of equations fl2"T| ) and ( P2"|) versus the separation r, figure 22 and 23. 

As a technical detail concerning the evaluation of the RKSH, we remark that L s = 
0.8 and that the dissipation fit is extracted from the range of separations [0, Ifu], with 
I f it = 0.2. According to our theoretical considerations, for r C L s we should recover the 
classical RKSH. However, the present value of 0.15 for S* implies that when r < L s we 
also have r ~ rj, hence no room is left for the establishment of the classical form of scaling 
below L,. 
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6.4 Extended self similarity 

Before closing the section, it is worth commenting on the issue of ESS scaling in shear 
flows. Actually, the use of equation (|20"D requires some caution in the present conditions, 
where a superposition of different laws may be present. Clearly, a blind fit through the 
entire available range can extract an effective slope. It is much safer, instead, to consider 
the local slopes in detail, to understand to what extent a power law is really observed, 
and, more interestingly, the range of scales over which a certain value for the exponent 
is found. As an example, in figures 24 and 25, the local slopes, p = 6 and p = 9, for 
the shear flow are contrasted with those of homogeneous isotropic turbulence. In both 
cases at small separation the dissipative scaling SV P oc r p is achieved confirming the good 
quality of the data. As we move towards inertial separations the local slopes of isotropic 
turbulence remain quite constant over almost one decade. For the shear flow instead a 
definite deviation from constant is found at approximately r/rj — 15, with the shear scale 
falling near by (L s /r) = 16). Similarly to what found in the lower part of the logarithmic 



region of the turbulent channel p6| , the local slopes do not display clear ESS scaling 
regions, presumably due to the competition of two different statistical trends and to the 
lack of sufficient scale separation to distinguish between the two. 

6.5 An estimate of the flatness 

A different estimate of the flatness is here proposed taking into account the similarity laws 
given by eq. fl2~T|) and (|22"D . Coming back to the issue of intermittency, let us consider 
again the flatness factor of the random variable 5V(r), 

< 5V 4 (r) > 

F(r) = "J 2 - 23 

< 5V 2 (r) > 2 



According to the two forms of RKSH, eqs. (|2l|) and (|22|) , the same flatness factor can be 



expressed in terms of moments of the dissipation field as 



< ef/ 3 > <e 2 r > 

Fb = 7 J/ 3 ^2 Fw = < e > 2 ' ^ 

< €r > ^ c r ^ 

respectively. The plots in figure 26 show the curves corresponding to the three equations 
(p3[) and (124]) . The values computed using the new RKSH, F w in eq. (|24]) , are in excellent 
agreement with those evaluated from the definition (|23|). It should be noted that the 
use of the wrong form of RKSH, corresponding to expression F b in eq. (|24]) , substantially 
underestimates the intermittency. 

The same analysis can be performed for homogeneous and isotropic turbulence. In 
the plots of figure 27, we show F, Ft, and F w for turbulence measured at the center of 
a jet far downstream (R\ ~ 800) [Q. As one can easily observe, the new form of refined 
similarity hypothesis fails in homogeneous isotropic conditions, as we should expect. 
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The ability of the new RKSH to predict the amount of intermittency in the velocity 
increments from the pdf of the dissipation field is, in our opinion, a strong check on the 
validity of the proposed form of scaling. 

7 Statistical properties of the dissipation field 

Given its role in both the new and the classical form of RKSH, the study of the dissipation 
field may give further insight in the behavior of shear turbulence. To this purpose figure 
28 reports the self-scaling of the dissipation field in logarithmic scale. The solid line with 
slope a = 3.1 shows clearly that for sufficiently large separations r, which correspond to 
small values of the two moments, the power law < > oc < > a is able to fit the data. 
The scaling is strikingly similar to that found in homogenous and isotropic turbulence, 
where the estimated exponent is 3. The self-scaling seems to hold for all the moments we 
have been able to check, and the values found for the shear flow are always in agreement 
with those of homogeneous and isotropic turbulence, see table |I|. This analysis suggests 
that relations of the form 

< 4 > oc <e 2 r > a " (25) 

with, more or less, fixed exponents are able to describe the statistics of the fluctuations 
of the dissipation field either in shear or isotropic turbulence. Hence, basic statistical 
properties of the dissipation field remain substantially unchanged in the two cases. This 
issue becomes even more evident when considering the second moment of the dissipation 
as a function of r. This quantity is reported in figure 29, where the triangles and the cir- 
cles, corresponding to isotropic and shear turbulence respectively, indicate a substantially 
identical behavior. All together, these findings suggest a universal form for the pdf of the 
dissipation field, p(e r /e) normalized by its mean value e which is obviously dependent on 
the flow field. 



8 Concluding remarks & comments 

The direct numerical simulation of a homogeneous shear flow in a confined box has been 
analyzed to address the issue of scaling laws and intermittency in the simplest conceivable 
environment where anisotropy is prevailing. The dynamics of the vortical structures 
observed in the DNS is characterized by a regeneration cycle associated to a substantial 
fluctuation in the level of turbulent kinetic energy. The regeneration mechanism has 
been discussed in detail, and our analysis confirms substantially to previous descriptions 
|T0|| which, however, where restricted to the initial phase of flow development. In fact, 



in complete agreement with [16| and |[L7|| , we confirm that the confinement of the flow 



to a periodic box generates a dynamics, which, although strongly characterized by the 
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regeneration cycle of the vortical structures, reaches statically stationary conditions. The 
characteristic time of the system seems to be quite insensitive to changes in the aspect 
ratio of the computational box, hence we feel that the present flow may be considered as 
a valid prototype for the study of shear flows. 

We find an increased intermittency, as evaluated by the flatness of the velocity in- 
crements, with respect to that observed in isotropic turbulence. In isotropic turbulence, 
the intermittency of the velocity increments is related to the intermittent behavior of the 
dissipation field via the Refined Kolmogorov Similarity Hypothesis. For the shear flow, 
instead, we observe a clear failure of the RKSH which confirms the conclusions previously 
reached from data in the buffer layer of a channel flow (DNS) and successively verified 
with hot-wire measurements in a flat plate boundary layer. 

To attempt a clarification of this point, we have considered the possible existence 
of a self-scaling behavior of the structure functions. Under this respect, the extended 
self-similarity expressing the p t/l -order structure function as a power law with exponent 
C(p) in terms of, say, the third order moment, is a well established feature of isotropic 
turbulence which is able to characterize the amount of intermittency in terms of the 
scaling exponents. For the homogeneous shear flow, as already observed in the lower 
part of the logarithmic region of turbulent channel flow, scaling exponents can be defined 
only for a limited region of scales (r < L 8 ) above which a competition between two 
different statistical behaviors makes it difficult to define any scaling at all. However, 
recent experimental results concerning the log-layer of a turbulent boundary layer at a 
higher Reynolds number, seem to indicate the existence of a double scaling regime [^7 . 



Nonetheless, we are presently able to grasp the nature of the intermittency in the shear 
flow, since the failure of the RKSH is accompanied by the establishment of a new form of 
scaling, as proposed in p] . Originally the idea behind the new scaling law was to account 
for the crucial role of the shear in originating a substantial production of turbulent kinetic 
energy via the Reynolds stress. In principle, this mechanism, directly associated to the 
presence of a mean shear, is independent of other details of the flow, such as the presence 
of nearby solid boundary or the specific dynamics of the vortical structures which sustain 
the turbulent fluctuations. This view is confirmed by the present results which show 
how the new form of RKSH holds not only for the buffer layer of wall bounded flows, 
as originally verified in JF], but also in the present homogeneous shear flow. Here the 
dynamics of the vortical structures, a part from a general resemblance, is in fact much 
different from that observed in the wall region and no wall is present to inhibit wall normal 
motions. Nonetheless the new scaling law is well established, confirming that its presence 
is the signature of a predominant effect of the shear in the proper range of scales. 

Beyond the direct observation of the new scaling law, the detailed analysis presented 
in the paper shows how the increased intermittency of the longitudinal velocity increments 
in shear dominated flow is, in fact, to be ascribed to the new scaling law. Actually, the 
statistical properties of the energy dissipation field are substantially similar to those we 
find in isotropic conditions. From our data, the scaling behavior of the energy dissipation 
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field can not be consistent both with the observed level of intermittency of the veloc- 
ity structure functions and with the existence of the classical RKSH. On the contrary, 
the new RKSH is able to consistently predict the level of intermittency in the velocity 
structure function in terms of the sole dissipation field. This confirms that the intermit- 
tency of velocity structure functions at shear dominated scales is largely affected by the 
instantaneous process of turbulent kinetic energy production. In other words, since the 
dissipation field is approximately unchanged, to compute the increased intermittency of 
the velocity in presence of shear we have to include explicitly the production mechanism 
in the scaling law. 

These comments suggest that the source of intermittency in the shear flow is in prin- 
ciple substantially different from that operating in isotropic turbulence. To speculate a 
little further on this point, let us recall that, a part from the aspect ratio of the box, 
the confined homogeneous shear flow is defined by two dimensionless parameters, namely 
S* and S*. These two parameters control the position of the shear scale L s relative to 
the length scale of the large eddies, l d , and to the Kolmogorov scale, 77, respectively. A 
large value of S* implies a substantial range of scales where the shear is crucial to the 
dynamics. On the other hand, a small value of S* corresponds to an extended range 
of scales where the dynamics is purely inertial and where the effect of shear is largely 
irrelevant. One should expect the classical scaling laws of homogeneous and isotropic 
turbulence to emerge clearly in the limit of S* small within the range L s » r » rj. 
A major contribution of the present paper is to confirm the existence of a new form of 
scaling law which emerges in the limit of S* large, in the range Id » r » L s . Ideally, 
when both S* and l/S* are large, one should expect a coexistence of the two laws in the 
two different scaling ranges. In this ideal conditions, the intermittent behavior of velocity 
increments is induced by the energy injection associated to the turbulent kinetic energy 
production mechanism and is originated in the shear dominated range, where the scaling 
is totally different from the classical RKSH. However, once the energy is transferred to 
the smaller scales where direct injection due to shear is irrelevant, the classical energy 
cascade mechanism takes over. The intermittency of the velocity increments at the shear 
dominated scales appears to smaller scales which belongs to the classical inertial range as 
large scale modulations. Along the classical line of reasoning, the large scale modulations, 
through the cascade, are reduced to the intermittent behavior of the dissipation field. This 
is tantamount to postulating a certain universality of the probability distribution of the 
dissipation field, once, in the spirit of K62 theory, universal values for the anomalous 
corrections to the exponents of the structure functions are assumed in the classical part 
of the inertial range. These concepts are consistent with our numerical experiments with 
the homogeneous shear flow, even though the scale separation is not sufficient to observe 
both scalings simultaneously. 

To conclude we like to comment about the different sets of data available to us in 
the present analysis. As figure ?? shows, a large shear dominated range, i.e. S* large, is 
not frequently met in practice. An exception is the buffer region of wall bounded flows. 
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There, almost the entire range of available scales falls above L s , which for this case is also 
very close to the Kolmogorov length, rj, since S* ~ 1. These are optimal conditions to 
observe the new RKSH as a pure law, without interference with its classical counterpart 
since the contiguous purely inertial range below L s does not exist at all. The logarithmic 
region, both from DNS and experiments, present L s ~ l d) i.e. S* ~ 1, and L s » rj. 
Here only the classical scaling can be expected. The conditions of the present confined 
homogeneous shear flows are somehow intermediate between the ones of the buffer and 
the log region respectively. Here L s falls almost exactly in the middle of the available 
range of scales, a more complex situation and by far the more interesting to discuss. If the 
shear dominated and the purely inertial range were large enough we might have observed 
the simultaneous presence of the two forms of scaling law. However this cannot be the 
case with the Reynolds number we can reach by DNS. In practice, what we observe in 
the present case can be described as a sort of superposition of the two regimes, with a 
dominating contribution from the new RKSH. 
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Table 1: Scaling exponents, a q , 
different, non-integer, values of q 
Horn. Iso. Horn. Shear 



of the moments < 
see eq. (|25|). 



e q r > with respect to < e 2 r >, for 



«2/3 
«4/3 
«5/3 
«6/3 
"7/3 
"8/3 
"9/3 



-0.11 
0.22 
0.55 
1.00 
1.56 
2.10 
3.00 



-0.11 
0.22 
0.55 
1.00 
1.57 
2.25 
3.10 
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